Bosonization of one dimensional fermions out of equilibrium 
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Bosonization technique for one-dimensional fermions out of equilibrium is developed in the frame- 
work of the Keldysh action formalism. We first demonstrate how this approach is implemented for 
free fermions and for the problem of non-equilibrium Fermi edge singularity. We then employ the 
technique to study an interacting quantum wire attached to two electrodes with arbitrary energy 
distributions. The non-equilibrium electron Green functions, which can be measured via tunnel- 
ing spectroscopy technique and carry the information about energy distribution, zero-bias anomaly, 
and dephasing, are expressed in terms of functional determinants of single-particle "counting" op- 
erators. The corresponding time-dependent scattering phase is found to be intrinsically related to 
"fractionalization" of electron- hole excitations in the tunneling process and at boundaries with leads. 
Results are generalized to the case of spinful particles as well to Green functions at different spatial 
points (relevant to the problem of dephasing in Luttinger liquid interferometers). For double-step 
distributions, the dephasing rates are oscillatory functions of the interaction strength. 

PACS numbers: 73.23.-b, 73.40.Gk, 73.50.Td 
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I. INTRODUCTION 

One-dimensional (ID) interacting fermionic systems 
show remarkable physical properties. The electron- 
electron interaction manifests itself in a particularly 
dramatic way in ID systems, inducing a strongly cor- 
related electronic state - Luttinger liquid (LL)i*2i2iii&. 
A paradigmatic experimental realization of quantum 
wires are carbon nanotubes 6 ; for a recent review see 
Ref.0 Further realizations encompass semiconductor—, 
metallic 9 and polymer nanowires 1 ^, as well as quantum 
Hall edges^ii 2 -. 

While equilibrium LL have been extensively explored, 
there is currently a growing interest in non-equilibrium 
phenomena on nanoscale and, in particular, in non- 
equilibrium properties of quantum wires. In a recent 
experiment^ 3 , the tunneling spectroscopy of a biased 
LL conductor has been performed (see also a related 
work on carbon nanotube quantum dots 14 ). A sim- 
ilar approach was used to study experimentally non- 
equilibrium quantum Hall edgesi^. Quite generally, the 
tunneling spectroscopy technique allows one to measure 
the non-equilibrium Green functions G< (r) . Analogous 
experiments^ have been carried out earlier in order to 
study energy distribution function and inelastic relax- 
ation processes in quasi-one-dimensional diffusive metal- 
lic samples. The interpretation of the results for a metal- 
lic sample is based on the Fermi liquid theory, and, in 
particular, on a kinetic equation for a quasi-particle dis- 
tribution function. In fact, even in that case, careful 
analysis requires taking into account non-equilibrium de- 
phasing processes^ which lead to additional broadening 
of the measured Fermi-edge structures in the tunneling 
current. In the case of strongly correlated, non- Fermi- 



liquid systems (such as LL) out of equilibrium, the situ- 
ation is much more complex. In this situation not only a 
quantitative theoretical analysis of , but even the very 
notions of quasiparticle energy distribution and dephas- 
ing, become highly non-trivial. The goal of the presented 
work is to construct a corresponding theory. To achieve 
this goal, we develop a formalism of non-equilibrium 
(Keldysh) bosonization. While we consider systems of 
ID interacting electrons in this work, we expect that it 
will be an important step in understanding the properties 
of a broader class of systems — non-equilibrium quantum 
fluids in low-dimensions. This includes, in particular, 
systems of cold atoms, with either fermionic or bosonic 
statistics. 

The structure of the present paper is as follows: 
In Sec.|TT]we discuss possible experimental realizations 
of a non-equilibrium LL. 

In Sec. IIIII we develop a bosonization technique for 
non-interacting electrons away from equilibrium. Work- 
ing within the Keldysh non-equilibrium formalism, we 
derive the action of the bosonized theory. While this 
action is quadratic at equilibrium (which is the essence 
of conventional bosonization), it now includes arbitrary 
powers in the bosonic fields. We demonstrate how this 
action can be used to express the Green function of 
non-interacting fermions in terms of a Fredholm func- 
tional determinant of a single-particle "counting" oper- 
ator (which is of Toeplitz type). We further discuss the 
relation between this problem and that of counting statis- 
tics. Specifically, our result is expressed in terms of the 
determinant at the value of the phase ( "counting field" ) 
A = 2ir. On the other hand, the counting statistics at this 
point is trivial, in view of charge quantization. We show 
that the difference between the determinants used for ex- 
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pressing the Green functions and those used for counting 
statistics results from different continuations (analytic vs. 
periodic) of the functional determinant beyond the non- 
analyticity point A = ir. 

In Sec. IIVI we apply our technique to the problem of 
Fermi edge singularity (FES) out of equilibrium. We 
show that non-equilibrium FES Green function is ex- 
pressed in terms of the same functional determinant but 
with a shifted value of the argument, A = 2(ir — So) , 
where <5o is the scattering phase on the core hole. Com- 
paring our results for this problem with those obtained 
earlier—, we establish useful identities between Fredholm 
determinants of counting operator at values of the count- 
ing field A differing by 2ir. 

In Sec|V] our formalism is extended to interacting 
fermions in a quantum wire. First, we analyze the prob- 
lem of tunneling spectroscopy of a non-equilibrium LL 
in the case of spinless fermions. We demonstrate that 
the non-equilibrium Green functions are expressed in 
terms of products of single-particle Fredholm determi- 
nants. The corresponding values of the counting fields are 
shown to be related to "fractionalization" of particle-hole 
excitations created during the tunneling process, as well 
as at the boundaries with non-interacting leads. Our re- 
sults for contain all information about single-particle 
properties of the system, including tunneling density of 
states, energy distribution, and dcphasing. We find, in 
particular, that the dephasing rate oscillates as a function 
of the interaction strength (LL parameter K), vanishing 
at certain values of K. At the end of the section we gener- 
alize the consideration to the case of spinful fermions, as 
well to Green functions at different spatial points (which 
is relevant to the problem of dephasing in LL interferom- 
eters). 

Section I VII includes a summary of our results as well 
as prospects for future work. 

Some of results of this work were presented in a Letter, 
Ref. GJ. 



II. NONEQUILIBRIUM LUTTINGER LIQUID: 
SETUPS 

In this section we specify the class of problems to be 
considered and discuss possible experimental setups. We 
assume that electrons with distributions functions n v (e) 
(rj = R, L labels right- and left-movers) are injected into 
a LL wire from two non-interacting electrodes. It is con- 
venient to model the electrodes as non-interacting ID 
systems, so that the whole structure is a wire with spa- 
tially dependent interaction that switches on near the 
points x — ±L/2, see Scc.|V]for detail. 

It is worth noting that we assume the absence of 
electron backscattering due to impurities inside the LL 
wire. When present in sufficient amount (so that one can 
speak about a disordered LL), such impurities strongly 
affect the electronic properties of a LL wire. Specifi- 
cally, they induce diffusive dynamics at sufficiently high 



temperature T and localization phenomena proliferat- 
ing with lowering T (Ref. [UHH), as well as inelas- 
tic processes^. We also neglect the nonlinearity of the 
electron dispersion whose influence on spectral and ki- 
netic properties of ID electrons was recently studied in 
Refs. [H |H. 

We discuss now possible experimental realizations of 
the problem. The simplest way to take the system out of 
equilibrium is to apply a voltage between two electrodes, 
so that the incoming distribution functions have different 
chemical potentials, fiL — fJ-R. = eV, but equal tempera- 
tures, Tr = Tl = T, see e.g. Ref. |2y. However, in the 
case of a LL this situation is almost identical to the equi- 
librium one, in view of the absence of electron backscat- 
tering. Indeed, the bosons remain at equilibrium, so that 
the usual bosonization technique (within Matsubara for- 
malism) can be applied. The only non-equilibrium effect 
will be a simple shift of the chemical potential of left- 
movers as compared to that of the right-movers. 

A generalization of this setup that does yield a non- 
trivially non-equilibrium LL is shown in Fig. [IJl. A long 
clean LL is adiabatically coupled to two electrodes with 
different potentials, fiL — ^r = eV and different tempera- 
tures Tl, Tr. (A particularly interesting situation arises 
when one of temperatures is much larger than the other, 
e.g., Tl — and Tr finite, so that non-equilibrium effects 
are most pronounced.) This model has been investigated 
in our previous works, Refs. I27ll28l . While showing gen- 
uinely non-equilibrium effects (in particular, energy re- 
distribution of electrons), this model, when treated in the 
framework of Keldysh bosonization formalism, is charac- 
terized by a Gaussian action. For this reason, we termed 
this setup "partially non-equilibrium" in Ref. (27|. We will 
verify in Sec|V]that the results of the present work (per- 
taining to full non-equilibrium) reduce to those obtained 
earlier (partial non-equilibrium) in the case when both 
hr and til are taken to be Fermi-Dirac functions. 

The focus of this work is generic non-equilibrium sit- 
uations, when at least one of the functions n v is not of 
the Fermi-Dirac form. Such situations naturally arise 
when electrons injected into a LL wire represent juxta- 
position of particles originating from reservoirs with dif- 
ferent chemical potentials and mixed by impurity scat- 
tering. Two possible realizations of such devices are 
shown in Fig. [TL,c. In the first case, Fig. Qi>, the mix- 
ture of left and right movers coming from reservoirs with 
fJ-L i= /J-r is caused by impurities which are located in 
the non-interacting part of the wires^. In the second 
setup, Fig.[lJ, the LL wire is attached to two thick metal- 
lic wires which are themselves biased. We assume that 
those electrodes are diffusive but sufficiently short, so 
that energy equilibration there can be neglected. As a 
result, a double-step energy distribution is formed in the 
electrodes^ and "injected" into the LL conductor. Such 
double-step distributions are of particular interest for our 
problem, as they are of the "maximally non-equilibrium" 
form. The existence of multiple Fermi edges in the dis- 
tribution functions "injected" from the electrodes ren- 
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FIG. 1: Schematic view of experimental setups for tunnel- 
ing spectroscopy of a LL out of equilibrium: (a) "partially 
non-equilibrium" setup, with distribution functions n v (e) of 
Fermi-Dirac form but with different temperatures; (b), (c) 
"fully non-equilibrium" setups characterized by double-step 
distribution functions n v (e) of electrons injected into the LL 
wire. 



ders the electron-electron scattering processe d ' which 
govern the non-equilibrium dephasing rate (and thus 
the broadening of tunneling spectroscopy characteristics) 
particularly important. 

The question of non-equilibrium dephasing induced by 
electron-electron scattering is particularly intriguing in 
the case of a ID system. First, energy relaxation is 
absent in a homogeneous LL system. Second, recent 
analysis of dephasing in the context of weak localiza- 
tion and Aharonov-Bohm oscillations has given qualita- 
tively different results: while the weak-localization de- 
phasing rate vanishes in the limit of vanishing disorder—, 
the Aharonov-Bohm dephasing rate is finite in a clean 
LL 22 ' 31 . In the case of a partially non-equilibrium setup 
the tunneling spectroscopy dephasing rate has a form 
similar to the equilibrium Aharonov-Bohm dephasing 
rat o 27 i 28 . As we show here, in the case of double-step dis- 
tributions dephasing acquires qualitatively distinct fea- 
tures; in particular, the dephasing rate becomes an oscil- 
latory function of the interaction strength. 

Having described the problems to be addressed, we 
turn to the corresponding formalism. It is instructive to 
develop it first for the case of non-interacting fermions 
and then "turn on" the interaction. 



III. FREE FERMIONS 

In this section we develop a bosonization formalism 
for the case of free fermions out of equilibrium. Specifi- 
cally, we consider non-interacting fermions with a given 



distribution function n(e) and derive the correspond- 
ing bosonic action. Using the latter, we calculate the 
fermionic Green function. Clearly, the Green function of 
non-interacting fermions is trivially obtained within the 
fermionic formalism. However, the results of this section 
are not just a complicated way to calculate a simple quan- 
tity. Rather, they will play a crucial role for developing 
the bosonic formalism for interacting systems studied in 
the remainder of the paper. 

A. Keldysh action: From fermions to bosons 

Bosonization has been proved to be a very effi- 
cient tool for tackling one dimensional problems at 
equilibrium^! 3 -!^, as it maps a system of interacting 
fermions (LL) onto that of non-interacting bosons. One 
can thus hope for similar advantages of this approach for 
non-equilibrium problems as well. The question though 
is whether the bosonization procedure can be generalized 
to systems out of equilibrium? As we show below, the 
answer is affirmative, yet substantial modifications are 
required. 

Quite generally, operator bosonization procedure con- 
sists of the following steps: (i) mapping between the 
Hilbert space of fermions and bosons; (ii) construc- 
tion of the bosonic Hamiltonian Hb representing the 
original fermionic Hamiltonian Hp in terms of bosonic 
(particle-hole) excitations, i.e. density fields; (iii) ex- 
pressing fermionic operators in the bosonic language; 
(iv) calculation of observables (Green functions) within 
the bosonized formalism by averaging with respect to 
the many body bosonic density matrix (ps)- Nei- 
ther the Hilbert space nor the operators (including the 
Hamiltonian) contain an information regarding a state 
of the many-body system. Therefore, the first three 
steps remain unchanged for a non-equilibrium situation. 
The major modifications occur in the step (iv). In- 
deed, at equilibrium the fermionic density matrix is ex- 
pressed through the corresponding Hamiltonian as pp = 
cxp(—Hf/T), implying that the same relation holds in 
the bosonized theory, ps — exp(— Hb/T), which makes 
averaging with respect to ps straightforward. Out of 
equilibrium this is not so anymore: a one-particle den- 
sity matrix corresponding to a non-equilibrium occupa- 
tion n(e) of fermionic states translates into a complicated 
density matrix of bosons, which does not allow the ap- 
plication of Wick theorem. This poses a major difficulty 
in bosonizing fermionic problems away from equilibrium 
and, as we see below, results in a non-gaussian action of 
the bosonized theory. 

To construct the effective bosonic theory, we start with 
the fermionic description. Within the LL model, the 
electron field is decoupled into a sum of left- and right- 
moving terms, 

1>(x, t) = i, R (x, ty^ x + M*, t)e~ VFX , (1) 

where pp is the Fermi momentum. The Hamiltonian of 
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the system reads 

H = -iv / dx (^dxipR - ip^dxipi 



(2) 



where v is the electron velocity. The bosonic representa- 
tion for fermionic operators has the formi^i 3 -^! 3 ^ 



(— 

\2irv 



1/2 



(3) 



where A is an ultraviolet cut-off. The bosonic fields § n (x) 
are related to the density of electrons (given by p v (x) = 
ipj ] (x)ip v (x) in the fermionic language) as 
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Pn\ x > = 2- d x<t>ri ■ 



(4) 



and obey the commutation relations 

[(/>r(x),(I>r(x')] = -[<t>L{x),<i>L{x')] = iTrsgn(x-x') . (5) 

We use the convention that in formulas rj should be un- 
derstood as r] — ±1 for right/left moving electrons. The 
bosonized Hamiltonian is expressed in terms of density 
fields in the following way: 



H = irv / dx (p R + p L ) 



(6) 



We turn now to the Lagrangian formalism. Since we 
deal with a non-equilibrium situation, the system is char- 
acterized by an action defined on the Keldysh contour 33 , 



S [tp]=Jdt J dx ^2 ^l id n^v 



(7) 



t)=R,L 



where ip, ip'(t, x) are fermionic fields, and 8r,l = dt ± 
vd x . To generate correlation functions, it is convenient 
to introduce a source term. 

S v {^] = j dt J da%(M)^ (MM, («,«)• (8) 

The field components on the upper branch and lower are 
denoted by + and — respectively. It is convenient to 
perform a rotation in Keldysh spaced 3 -, thus decomposing 
fields into classical and quantum components (the latter 
being denoted by a bar), 



V r , ) V n = (V + , t ,±V- jV yV2, (9) 

PviPv = (P+,v ± P-,v)/^2> ( 10 ) 
i'n^rj = (ip+v ± i>- lV )/V2, (11) 

i^ = (<,T<,)/V2- (12) 

In these notations, the density correlation functions are 
encoded in the generating function 



Zr,[Vr,, V v ] = (expliVvPr, + iVr,Pr,})So 



(13) 



The calculation of the partition function can be per- 
formed in cither the fermionic or the bosonic description. 
In the fermionic language it can be readily done by eval- 
uating a Gaussian integral over the Grassman variables, 



Z V [V,V} 



det[l + G n o(cro^ + <7itOA/2], (14) 



where ctq and a\ are the unit matrix and the first Pauli 
matrix in the Keldysh space, and G v o is the Keldysh 
Green function of free chiral fermions, which has the fol- 
lowing matrix structure: 



G 



77O 







/To, 



(15) 



Here G^ , GJJ and G^ are advanced, retarded and 
Keldysh components, 



T n0 



(e,p) = l/(e — r\vp ± iO) . 



(16) 



G* (e,p) = [1 - 2n v (e)][G r v0 (e,p) - G^(e,p)] .(17) 

We expand now the generating functional (I14j) in pow- 
ers of the source fields V n , V v . For higher-dimensional 
systems this would generate all terms of the type V™V™ . 
In ID the situation is different. Specifically, in an equi- 
librium ID system only terms up to second order (V^V^ 
and V^ ) are generated, which forms the basis of conven- 
tional bosonization. Out of equilibrium, this is not true 
anymore: terms of higher orders are generated as well, 
and the theory becomes non-gaussian. What is crucial, 
however, is that all higher-order terms are of the type 
V v , i.e. they do not depend on V v . We will prove this 
statement in Sec. MI Bl and ITTT"Cl below. 

The generating functional has thus the structure 



Z n % V] = exp -iVnU% + J2 l ~K S ^v > (18) 



where S n ,r] is the n-th order irreducible vertex function, 

*^n,?7 (^1 j ^1 j ■ ■ ■ 1 %n j in) 

= -i n ^2 r TrKG 1l0 (xi,t 1 ]x l21 ti 2 )^j= 



X GrjO {Xi 2 i ti 2 i x i 3 , ti 3 ) X 
X GriQ(Xi n ,ti n \X\,t{) 



V2' 



(19) 



The multiplication in Eq. (|18[) and analogous formulas 
below should be understood in the matrix sense with re- 
spect to the coordinates, 



VrffiVr, = J [dxWdt^X^h) (20) 

x n^ l {x 1 ,t 1 ]x 2 ,t 2 )V ri (x 2 ,t 2 ) , 

= J [dx] [dt] V v (xx, h) . . . V v (x n , t n ) 

x S n ^[x\^ t\\ . . . ; x n , t n j , (21) 



where J" [da:] [dt] implies integration over all spatial and 
time coordinates. The summation in Eq. (|19p goes over 
(n— 1)! permutations is, . . . , i n } of the set of indices 
{2, . . . , n} (labeling the space-time coordinates), and Ttk 
denotes the trace over Keldysh indices. Clearly, after 
integration with V v fields in Eq. (fT8|) all the (n — 1)! 
terms of the sum in Eq. (|19|) yield equal contributions, 
so that the total combinatorial factor is (n— l)!/n! = 1/n, 
as should be in the expansion of the logarithm. We 
have chosen to define the vertex function in the sym- 
metrized form (|19p [and to introduce the corresponding 
factor \/{n— 1)! in Eq. (|18[) ], since S n ^{x\, t\; . . . ; x n , t n ) 
are then equal to irreducible density correlation functions 
((/?(xi,ti) . . . p(x n ,t n ))) . 

The quadratic part of the generating functional (fT8|) is 
determined by the polarization operator of fermions, 



o n° 



K 



(22) 



with the retarded, advanced, and Keldysh components 
given by 



n r,a 



1 



rjq 



(23) 



27r 7yvg — u> =F zO 

n£( w , <?) = [ie>, 9 ) - il>, 9 )]b» . (24) 

Here the function 
1 



B,(w) = 



den v (e) [2 — n v (e — to) — n v (e + u>)], (25) 



governs the distribution function N v (co) of electron- 
hole excitations moving with velocity v in direction rj, 
B v (ui) = 1 + 2N v (ui). At equilibrium, 

B„(w) = S cq (cj) = 1 + 2iV cq H = coth (^) , (26) 

where N eq (to) is the Bose distribution. By construction, 
the second order density-correlation function S v . n =2 in 
Eq. (fig)) is equal to the Keldysh component ITr of the 
polarization operator (times —i). 

In order to bosonize the theory, we should find a 
bosonic counterpart of the action Sq^ that reproduces 
the generating functional (fTS|) . According to Eq. (fT3")) . 
we have 

exp (i5 0r) [^,^])= /" VV v V%Z n \y v ,V v ]e^ v *to-™'*»' . 

(27) 

Substituting Eq. (TTSJ) into Eq. l[2"7]). we obtain the 
bosonized action 

"Sor, [Pt, , Pr,] = P,,n^ 1 p I( - i In Z ?) [x] • (28) 

Here Z v [x] = Z v [x = 0,x] is a partition function l|18p of 
free fermions, 



(29) 




FIG. 2: Vacuum loops for free fermions in an external field 
V. At equilibrium only S2 is non-zero, according to the 
Dzyaloshinskii-Larkin theorem. Away from equilibrium, all 
vertices appear. For details see Sec. IIII Bl 



subject to the external quantum field 



(30) 



The combined action of left- and right- moving electrons is 
simply given by a sum of the corresponding chiral actions, 



So[p,p] = ^2S 0v [p v , 



Pnl 



(31) 



Thus we have described a system of non-equilibrium 
free fermions by a bosonic theory, Eq. ([3"Tj) . In this ap- 
proach information on the non-equilibrium state of the 
system is encoded in the vertices (S nv ), schematically 
depicted in Figf2] In Sec. IIII Bl we discuss the status and 
implications of the Dzyaloshinskii-Larkin theorem con- 
cerning these vertices. 



B. Dzyaloshinskii-Larkin theorem 

The appearance of higher-order (n > 2) fcrmionic ver- 
tices may seem to contradict the Dzyaloshinskii-Larkin 
theorem^. The latter states that diagrams containing 
closed loops with more than two fermionic lines vanish. 
Although the theorem was formulated for the equilibrium 
case, its proof, given in Ref. [34, ostensibly relies solely 
on particle conservation. Since the latter remains valid 
out of equilibrium, one might expect the theorem to hold 
under non-equilibrium conditions as well. To understand 
why Dzyaloshinskii-Larkin theorem is in fact restricted to 
the equilibrium case only, and what its implications for 
a non-equilibrium situation are, we carefully re-examine 
the arguments of Ref. [34]. 

One starts with the continuity equation for the chiral 
current and density operators, 



n=2 



up v - min = ■ 



(32) 
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Since within the LL model these operators are related to 
each other through j v — rjvp r)1 the continuity equation 
can be rewritten in terms of the density field only, 

(uj - ijvq)p v = . (33) 

As a consequence, correlation functions of densities sat- 
isfy 

(u>i - r)vq i )(p v (u>i,q 1 )p v (u} 2 ,q2) ■ ■ .p n (oJ n ,q n )) =0 (34) 

for any i — I, ... ,71. Therefore, the irreducible density 
correlation functions S nri (u>x, qx; u>2, q%\ ■ ■ ■ w n) q n ) with 
n > 2 should be zero everywhere, except possibly for 
the mass shell with respect to all arguments^, 

= 5{u)i - rjvq^Si^ - r)vq 2 ) . ■ . S(u> n - f]vq n ) 

xSr,(LJ-i,UJ2, ■ ■ ■ ,U n )S(uJl + ■■■ +W re ). (35) 

In the case n = 2 the argument is not applicable in view 
of the Schwinger anomaly, yielding the first term in the 
exponent on the r.h.s. of Eq. (|18[). When translated 
into the coordinate-time space, the mass-shell condition 
(|35[> implies that the correlation function depends in fact 
only on the world line to which each of the points (U, Xi) 
belongs but not on the position of this point on the line: 

S n ri(ti,x 1 ;...;t n ,x n ) = {{p v (h,Xi)...p TI (t n ,x n ))) 
= «^(0,£i)...p„(0, £„))>, 

(36) 

where £j = x — rjvti. In the Keldysh formalism language, 
the only non-zero irreducible density correlation function 
in any order n > 2 arises when one considers the corre- 
lator with all n fields being the classical components p. 
(This follows from the fact that the operators p commute 
to a c-number.) These correlation functions are the noise 
cumulants in the system. 

The behavior of the correlation functions S nri on the 
"light cone" ((35)) can not be determined from particle 
conservation law and requires an additional calculation. 
While at equilibrium all S nr] with n > 2 do vanish (which 
reconciles our theory with the Larkin-Dzyaloshinskii the- 
orem), out of equilibrium they are in general non-zero. 
We consider this general situation in Sec. IIII Cl where we 
show that the bosonized action can be presented in a 
compact form of a functional determinant. 

C. Bosonized action as functional determinant 

As we have shown, the bosonic action, Eq. (j3"Tj) . 
is expressed through the partition function Z[V, V] of 
free fermions in an external field V(x,t) defined on the 
Keldysh contour. In one dimension the partition func- 
tion can be cast in a relatively simple form. To achieve 
this, we first present the partition function 

Z[V,V] =tr{p F S c }. (37) 



Here S c is an evolution operator along Keldysh contour, 
Z[V,V] - Urn tl { pFe -mv + (-t)]At e -iH[v + (-t+At)]At 

t—*oo 

X ... X e -iH[V+(t)]At e iB[V-[t)]At e iH[V-(t-At)]At 

x ... x e mv-(-t)]At h (38) 

and the trace is taken over the many-body fermionic Fock 
space. Equation (|38[) can be further simplified by means 
of the following identity^ 

tr{e Hl e H2 ... e" N } = det (1 + e hl e h2 . . . e' lN ) . (39) 

Here h„ is a matrix in the single particle Hilbert space, 
and 

X /';;'";«, (40) 

is the corresponding operator quadratic in fermionic cre- 
ation/annihilation operators (a\a). The trace in the 
l.h.s. of Eq. (f3"9|) is taken in the many-body Fock space, 
while the determinant on the r.h.s. is taken in the single- 
particle space. 

Applying Eq. (|39[) in the continuum limit, we express 
the partition function in the following form 

Zr,[V v , V v ] = det[l -n n + r^U^U-a] ■ (41) 

Here 

U+^(t) = Tcxp ^— i J dth + ,^j , 

UZ^it) = Tcxp (i dth-,^j (42) 

are evolution operators that correspond to the single- 
particle Hamiltonians 

Q 

h- iV = -ir}v^ + V-{x,t). (43) 

Thus the many-body problem of summing all vacuum 
loops has been reduced to a calculation of a functional 
determinant of an operator in a single-particle Hilbert 
space. To simplify it further, we analyze the properties 
of the evolution operator U in one dimension. Its action 
on a wave function ip(x) can be described as 

ip(x, t) = Texp (-i dth^j i/j(x, 0) , (44) 

where ip(x,0) = ip(x). One can easily show that the 
resulting wave function %p(x,t) satisfies the Schrodinger 
equation 

d 

i—J>{x,t) = h+ij}{x,t) . (45) 
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Solving Eq. ([45l) explicitly one finds 

i){x, t) = i(j(x - nvt, 0)e- 1 Jo dTV + (x +v v(T-t),T) ^ 

Therefore, the action on a wave function of the evolution 
operator forward and backward in time results in the 
phase factor 

(EC 1 17+ V) (ar) = lp{x)e-^o dr{V+-V-){ x +vr,r) ( 4? ) 

Consequently, the partition function of the ID 
fermions can be cast as^ 

Z n [V v: %] = e^ n ^A ri [S ri (t)} , (48) 

where we introduced a determinant 



A,, [<*,,(*)] = det[l + (e~ 



l)n v ) 



(49) 



6 n (t) =V2 dTV v (r)v(T + 1), t) (50) 



and 



is the scattering phase accumulated by an electron mov- 
ing along a "light-cone" trajectory. Thus, according to 
Eq. (|48[) the problem of summing up the vacuum loops 
is reduced to evaluation of a one-dimensional functional 
determinant (|49|) . 

The determinant (|4"9"1) is defined by the function 8 v (t) in 
the time space and n v (e) in the energy space, with e and t 
understood as canonically conjugate variables. It belongs 
to the class of Frcdholm determinants. For a specific case 
(that will be particularly important for us below) when 
S v (t) is different from zero within a limited interval of 
time only, the determinant acquires the Toeplitz form. 
Such determinants have been considered in the context 
of counting statistic o 38 i 39 ; see a more detailed discussion 
in Sec. IIIID1 It is also worth mentioning that there is 
a vast literature on the connection of Fredholm determi- 
nants to quantum integrable models, classical integrable 
differential equations (with soliton solutions), and free- 
fermion problems; we refer the reader to Refsi 40 i 41 ' 42 i 43 
and references therein. 

At equilibrium the Taylor expansion of In Z in S ter- 
minates at the second order (S n — for n > 2), in 
agreement with Dzyaloshinskii-Larkin theorem, Ref. [34]. 
In that case the action (f3"Tj) is quadratic, reproducing 
the standard LL model. Away from thermal equilib- 
rium, high-order density correlations are finite 3 ^,. For 
this reason, we obtain a non-Gaussian bosonized theory, 
despite the fact that the Hamiltonian © is quadratic. 
The higher-order terms S n with n > 2 appear in the 
bosonic action due to a non-diagonal structure of the 
density matrix in the bosonic Fock space, which leads to 
a breakdown of Wick theorem for the bosonic fields. 



D. Green functions 

We have thus shown that non-interacting fermions can 
be equivalently described by the bosonic theory with the 



action given by Eqs. ([31]), (EH), We apply now this 

formalism to calculate the free-fermion Green functions 
(GFs), 

G<(x 1 ,ti;x 2 ,t 2 ) = i(ipl(x2,t2)ip v (xi,t 1 )), 
G>(x 1 ,t 1 ;x 2 ,t 2 ) = -i(ipr l (xi,ti)ifjl l (x 2 ,t 2 )} . (51) 

At equilibrium these GF's are related to the advanced 
and retarded GFs via 



G>(x, e) = [G r v (x, e) - G*(x, e)](l - n v (e)) , 
G<(x,e) = -[G r Jx,e)-G a Jx,e)]n ri (e). 



(52) 



For free fermions, Eq. (|52p is valid for an arbitrary dis- 
tribution function n v determining the filling of single- 
particle states. 

Due to Galilean invariance the GFs depend only on 
t v = t\ — t 2 — r/(xi — x 2 )/v, so we may set x\ = x 2 = 
vt 2 = in the argument of GF. Using Eqs. ([3]), ([5T]) . we 
obtain 



iA 
2ttv 



(T K e 



.(0,T„) e -i4v + (0,0)\ 



(53) 



and a similar result for the function Gq v . 



At thermal 



equilibrium Gq v can be readily calculated. A standard 
calculation (presented for completeness in Appendix IA"]) 
yields 



G}{r n ) = T 



2v sinh ttTt^ 1 ± iAr^ 



(54) 



Away from equilibrium the calculation of GF's, rather 
simple within a fermionic framework, turns out to be 
quite complicated within a bosonic one. Nevertheless, 
this effort pays off, since the bosonic formalism will later 
allow us to extend the analysis to the interacting case. 

Within the bosonic description the GF can be repre- 
sented as a functional integral over the density fields. 
Since calculations of Gq v and Gq ^ are quite similar to 
each other we focus here on 



iA 
2ttv 



VpVpe lS[pti] 

,)-0(O,O)-0(O,r, 



-0(0,0)] 



(55) 



In a generic non-equilibrium situation the bosonic ac- 
tion, Eq. ((31]) , contains terms of all orders with no small 
parameter; the idea to proceed analytically in a con- 
trolled manner may seem hopeless. This, however, is 
not the case: non-equilibrium bosonization is an efficient 
framework in which the functional integration can be per- 
formed exactly. Indeed, Z v in Eq. |28|) depends only on 
the quantum component p, so that the action, Eq. (|31|) . 
is linear with respect to the classical component p of the 
density field. Hence the integration with respect to p can 
be performed exactly 



iA f 

Go. v {t) = -— I VpZ r) [xr ) }5{dtp + 'nvd x p- j) 



x 



[0(O,r)+0(O,O)] 



(56) 
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where the source term is 

j(x,t)=6(x)[6(t-T)-6(t)]/V2. (57) 

Resolving the ^-function, we obtain an equation that de- 
termines the quantum component of the density field, 



dtPn + Wd x pr, = j(x, t) . 



(58) 



According to the structure of the first term in the ac- 
tion (|28|) . we should look for the advanced solution of 
Eq. (|58p which is zero at times larger than those at which 
the source j(x, t) acts. In other words, in the asymptotic 
regions \x\ > L/2 the solution p(x,t) should contain in- 
coming waves only. Solving Eq. ([55| with this asymptotic 
conditions, we find the quantum density component 



p v (x,t) 



9(—T]X) 

~7T~ 



S(x — r/vt) — S(x — rjv{t — r)) 



(59) 



To find the Green function, we need to evaluate the fac- 
tors multiplying the delta- function in Eq. (|56p . subjected 
to the 8- function constraint. The most non-trivial fac- 
tor (which carries the information about the distribution 
function) is Z^x-q], where Xn is related to p v via Eq. (|30[) . 
According to Eq. (|4"5)) . Z v [x v ] is expressed as a functional 
determinant of the form (l49l) . We thus obtain 



1 



1 



2ttv t =p if A 



(60) 



Here we have denoted by the determinant normal- 
ized to its value for zero-temperature equilibrium dis- 
tribution, see Appendix A. It is convenient to use this 
definition since the determinant A,, requires in fact an 
ultraviolet regularization. On the other hand, the nor- 
malized determinant A v (which is equal to unity for the 
equilibrium, T — case) is uniquely defined. The prefac- 
tor in Eq. (f6"0"| that does not depend on the distribution 
function is immediately determined from the equilibrium 
result. 

According to Eqs. (|4"51) , (|5U|) . the mass-shell nature of 
S vn implies that Z^lx^] depends only on the world- line 
integral 



5 v (t) = V2 / dt x v (vvt, i-t). 



(61) 



Using Eq. ([30|) . we find an explicit solution for the 
"counting field" Xr], 



X v (x,t) = 2ir[vp n (x,t) + 77 / dxp n (x,t)}. (62) 

Jo 

Next we calculate the value of S(t) for our non-interacting 
problem. Substitution of Eq. ([62)) into Eq. (|6"Tj) allows us 
to cast the result for the phases S^it) into the following 
form: 

drj^t) — ~27tv / 2 77 lim / dxp v (x,t). (63) 

{— >— 00 Jo 



For the free-fermion problem the phase 8 n (t) = 
Xu! T (t, 0) where 

w T (t,t) = 6(t - 1) - 6{t - t - t) (64) 

is a "window function" and A = 2tt. Thus, Z v [xrl\ — 
A J)T (27r), where A^ T (A) is the determinant (|4^1) (normal- 
ized to its T — value) for a rectangular pulse. 



^0,77 ( r ) 



1 A^(2tt) 
2~kv t =p i/A 



(65) 



Determinants of the type (|4"9")) have appeared in a the- 
ory of counting statistic a 38 ' 39 . Specifically, the generating 
function of current fluctuations k(X) — X)^L-oo emX Pn 
(where p n is the probability of n electrons being trans- 
ferred through the system in a given time window r) 
has the same structure as A^ T (A). Taylor expansion of 
lnK(A) around A = defines cumulants of current fluc- 
tuations. 

According to its definition, k(A) is 27r-periodic, which 
is a manifestation of charge quantization that should 
show up in measurements of the transfered electric 
chargo 37 i 38 i 39 i 44 i 45 i 46 i 47 . Thus, k(2tt) = 1 is trivial. On 
the other hand, we have found that the free electron GF 
is determined by the non-trivial value of the functional 
determinant exactly at A = 2n. A resolution of this ap- 
parent paradox is as follows: the determinant A T?T (A) 
should be understood as an analytic function of A in- 
creasing from to 2ir. On the other hand, k(A) is non- 
analytic at the branching points A = ±7r, ±3ir, .... To 
demonstrate this, it is instructive to consider the equi- 
librium case that is treated in Appendix fAJ Then the 
expansion of In A^ T (A) in A is restricted to the A 2 term 
(since RPA is exact). It is easy to check that the A = 2n 
point on this parabolic dependence correctly reproduces 
the fermion GF via Eqs. ([101), (US])- As to the count- 
ing statistics ln/c(A), it is quadratic only in the interval 
[— 7r, 7r] and is periodically continued beyond this interval, 
see Fig. H 

The difference in the analytical properties of k(A) and 
A^ T (A) becomes especially transparent if one studies the 
semiclassical (long-r) limit, 



In A^ T (A) = — / del \n[l+{e- tX -l)n v (e)}+t\6(-e) 

(66) 

For small positive A the singularity of the integrand clos- 
est to the real axis is located at e = z(tt — A)T, i.e. near 
e = iirT. As A increases, the singularity moves towards 
the real axis, crosses it at A = ir and finally approaches 
e = —iirT as A — > 2ir (see inset of Fig. [3]). The inte- 
gral for ln«(A) is taken along the real axis, resulting in 
non-analyticity at A = ir and in zero value at A = 2ir. 
On the other hand, the contour of energy integration for 
lnA r)T (A) with A > 7r is deformed in the complex plane 
to preserve analyticity, as shown in Fig. [3] Specifically, 
the contour consists of the integration along the real axis 



a part along the branch cut on the imaginary axis. The 
integration along the real axis yields 



de 



In 



e T 



-iX 



iX0(-e) 



TX 2 , s 
= (67) 



where N = [X/2tt], X = X + 2irN. The integration along 
the branch cut of the logarithm yields -(T/2)[(27riV) 2 + 
47riVA], resulting in the long-r asymptotics 



lnA J)r = -tTA 2 /4tt. 



(68) 



Substituting this in Eq. I|65p , we correctly reproduce the 
long-time asymptotics of the Green function Gq at equi- 
librium, Eq. ([54]) . 

Let us now turn to the non-equilibrium situation and 
consider the double step function 



n v (e) = a n n Q (e_) + (1 - a n )n Q (e + ) . 



(69) 



where no(e) = 6(— e) is the zero-temperature Fermi-Dirac 
distribution function, e± = e — \i ± V/2, and < < 1. 
The value of [i is fixed by demanding that the total num- 
ber of electrons is the same as for the equilibrium distri- 
bution iM)(e) (which we use for normalization), yielding 
e_ = e — (1 — a)eV, e + = e + aeV. The distribution 
function in the time domain 



n v( T ) 



de 
to* 



-ier+Oe 



n v (e) 



(70) 



can be straightforwardly calculated, and is given by a 
sum of oscillating terms 



n v (r) = (1 - a,)e^ eyT no(T) 



where 



no(r) 



1 



2n t + iO 



(71) 



(72) 



is the T = Fermi-Dirac distribution function in time 
representation. 

On the other hand, we can find the time dependence 
of the fermionic distribution function by using our non- 
equilibrium bosonization approach, leading to the iden- 
tity (|65p. In the long-time limit we need to evaluate the 
integral (|66|) . yielding 



hiA, )r (A) 



eVr 
~2^~ 



( In (l - a, ; + a n e lA ) + aiX\ 



Analytically continuing in A we get 

1 , a v > 1/2 
a v < 1/2, 



lnA r;T (27r) ~ ieVr 



(73) 



(74) 



which reproduces the long-time time limit of the Green 
function of free fcrmions with the distribution function 



2 





17tT 




i7tT 




-In 



FIG. 3: Analytic A,, T (A) vs. periodic k(A) continuation of 
the functional determinant. The value of A VT at A = 2tt 
determines the free-electron GF, while In k(2tt) = in view 
of charge quantization. As an example, the equilibrium case 
is shown. Inset: contour of integration for the quasiclassical 
limit, Eq. (|66[) . of A,, T (A) is deformed, since a singularity of 
the integrand crosses the real axis at A = n. 



(fTTj) . We have just demonstrated how the identity 
works for a double-step non-equilibrium distribution. 

Equation (|65[) is a remarkable identity, as it connects 
two seemingly unrelated objects: the distribution func- 
tion of free fermions and a Fredholm determinant of the 
counting operator. The value of A = 2-7T appearing in the 
bosonic representation of the free-fermion GF Go,r ) (r) has 
a clear physical meaning: a fermion is a 27r-soliton in the 
bosonic formalism. 



IV. FERMI EDGE SINGULARITY 

A natural question to ask is whether values of A T?T (A) 
away from A — 2ix are physically important. To see that 
this is indeed the case, consider the Fermi edge singular- 
ity (FES) problem. In this problem, an electron excited 
into the conduction band, leaves behind a localized hole, 
resulting in an s-wave scattering phase shift, of the 
conducting electrons^. In the mesoscopic context^, the 
FES manifests itself in resonant tunneling experiments^. 
On a formal level it is described by the following Hamil- 
tonian 

ia\ak'bbK (75) 

While in the FES problem there is no interaction be- 
tween electrons in the conducting band, it has many fea- 
tures characteristic of genuine many-body physics. His- 
torically, the FES problem was first solved by an exact 
summation of an infinite diagrammatic series^. Despite 
the fact that conventional experimental realizations of 
FES are three-dimensional, the problem can be reduced 
(due to the local character of the interaction with the core 
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hole) to that of one-dimensional chiral fermions. For this 
reason, bosonization technique can be effectively applied, 
leading to an alternative and very elegant solution^ 1 -. 

Away from equilibrium, the FES has been addressed in 
Ref. [l8| where the canonical (fermionic) FES theory was 
combined with the scattering matrix approach. Below, 
we apply the non-equilibrium bosonization technique to 
the same problem. 

As mentioned above, the FES problem is effectively 
described by chiral ID electrons interacting with a core 
hole that is instantly "switched on". As was shown in 51 , 
taking into account the core hole in the bosonization ap- 
proach amounts to replacement of e 4 * by e l 0-- s o/T)4> j n 
the boson representation of the fermionic operator. Us- 
ing Eqs. JIT]), one gets 



G>(r) 



iA 
2nv 



/ Tife i(i-^>-(o,T) e -i(i-ia)0 + (o,oh 



(76) 



and similarly for the function G < . Within our non- 
equilibrium formalism, this implies a replacement j — > 
(1 — So/ir)j in Eq. (|58[) . Performing the derivation as in 
the free-fermion case, we thus obtain the non-equilibrium 
FES GF for electrons with an arbitrary distribution n(e), 

G*(t) = T«AA r (2 7 r-2£ )/27™(l±iAT)( 1 -' 5o / 7r)2 . (77) 

At equilibrium Eq. ([77)1 can be further simplified (see 
Appendix A), 



G*(r) 



(- 



_i0 )2 



=rA 



V sinh ttTt I 



reproducing the known result a 48 ' 51 . 

For a double-step distribution, Eq. 
limit is obtained as 



A r (27r-2<5 ) ~ e - T/2T f {2So) 
where the dephasing rate tT 1 is given by 



(78) 
the long-time 

(79) 



with the closed-loop part given by 
e c = A T (-2tf ) 



(82) 



This representation of the Green function is similar to the 
functional bosonization approac h 52 i 53 ' 54 i 55 , that employs 
both fermionic and bosonic variables. While functional 
and full bosonization approaches yield equivalent results, 
this equivalence is highly non-trivial. Indeed, comparing 
Eq. ||5TJ) with (J77J) and employing Eq. ([55]). we establish 
the identity 



=F 



iA 



(1-5o/tt) 2 



2irv \l±iArJ n 07 

= L*(t)A t (-25 ) 



(83) 



relating the functional determinants A r (27r — 2Sq) and 
A r (— 2 Jo) through the single-particle Green function 
L(r). 

Since is diagonal in energy space, while S v is di- 
agonal in time space, they do not commute, making the 
determinant non-trivial. It is worth noting that the func- 
tional determinants A r (A) for |A| < ir have been effi- 
ciently studied by numerical mean a 56 i 57 . The identity 
([83]) can be useful for the numerical evaluation of A r (A) 
at larger values of A. 



V. INTERACTING ELECTRONS 

So far we have been dealing with non-interacting elec- 
trons. Now we focus on the main subject of this work: 
bosonization of interacting fermions, both for spinless 
and for spinful cases. We begin by showing in Sec. IV Al 
how the interaction can be incorporated into the non- 
equilibrium bosonization scheme developed above. 



A. Keldysh action 



"(A) 



•— In ( 1 - 4a(l - a) sin 2 - 
2tt \ 2 



(80) 



In the energy representation t^ 1 determines the broad- 
ening of the split FES singularities. The same result for 
the broadening of FES has been obtained by Abanin and 
Levitov in Ref. [l8| within the fermionic framework. It 
is instructive to compare their result with our analysis. 
In the bosonization technique we have expressed the GF 
of the FES problem in terms of a functional determinant 
([77)) . On the other hand, within the fermionic approach 18 
the GF splits into a product of an open line L(t) (i.e. sin- 
gle particle Green function of fermions in the presence of 
external time dependent field) and closed loop e c (i.e. 
vacuum loops of fermions in an external field) , 



G*(t) =L^{r)e c 



(81) 



For the problem of spinless interacting fermions the 
Hamiltonian reads 

H = H +H CC1 (84) 
where Hq is given by Eq. @ and 

H cc = ^ J dxg(x) {p L {x) + p R {x)f , (85) 

where g(x) is a spatially dependent interaction strength. 
To model the coupling with non-interacting leads, we 
will assume that g(x) is constant within the interacting 
part of the wire and "switches off" near the end points, 
x = ±L/2, see Fig. [4] This way of modeling leads was 
introduced in Refs. 58 59,60 to study the conductance of 
a LL wire; it was also exploited in Refs. H2-62 to analyze 
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FIG. 4: Schematic view of a LL conductor connected to leads 
with two different incoming fermionic distributions. The LL 
interaction parameter K(x) is also shown; the dashed line 
corresponds to its sharp variation at the boundaries. 



the shot noise. In the Lagrangian formulation, Eqs. (|84|) . 
correspond to the action 



s[i>] = SoM + Sc#], 

S[ip] — j ' dt J dx^^ ^jjid^rf 



/0) 2 



p\x,t) 



where p{x) = Pr[x) + Pl(x). Decoupling the interaction 
term via a bosonic field tp by means of a Hubbard-Stra- 
tonovich transformation, we obtain the action 

S[ip,tp] = J^dt J dx ^2 ( id v- vWv 



■q=R,L 



-ipg 1 (x)ip 



(86) 



The theory of fermions in an arbitrary field tp(x, t) (on the 
Keldysh contour) can be bosonized using the results of 
Sec. Mil Introducing, as before, notations with (without) 
bar for the quantum (classical) components, we obtain 
the action 



S[<P,<P,P,P] = So[p,p] + S ee [p,<p] 



(87) 



where Sq [p, p] is the bosonized action of non-equilibrium 
free fermions, Eq. (f3"Tj) and 



See[p, f] 



dtdx{tpp + tpp + <pg (x)(f] 



Integrating out the auxiliary Hubbard-Stratonovich field 
tp, we derive a theory written solely in terms of density 
fields, 



S[p, p] = S [p, p]- J dtdxg(. 



x)pp. 



(89) 



Equation (|89|) constitutes a bosonic description for inter- 
acting electrons out of equilibrium. 



B. Tunneling spectroscopy of interacting fermions, 
spinless case 

We are now prepared to address the problem formu- 
lated in the beginning of the paper: an interacting quan- 
tum wire out of equilibrium, Figf3J We will first calculate 



the GFs at coinciding spatial points, which corresponds 
to tunneling spectroscopy measurements. In Sec. IV Cl we 
will generalize this analysis to GFs at different spatial 
points which are, in particular, relevant to experiments 
on LL interferometers. 



1. Tunneling into the interacting part of the wire 

We consider G\{t) for the tunneling point (x = 0) 
located inside the interacting part of the wire (region II 
in Fig. |4]); generalization to tunneling into one of non- 
interacting leads (regions I and III in Fig. Q| is straight- 
forward and will be presented in Sec. IV B 31 

Proceeding in the same way as for the non-interacting 
case, we come to a representation of the GF in the form 
of an integral over the density fields p and /?, Eq. (|55|) . 
The only difference as compared to the non-interacting 
case is that the bosonic action (j89[) now contains also the 
second term induced by the interaction. Since this term 
is linear in the classical component p v , we can perform 
the integration over it in the same way as we did in the 
non-interacting case. As a result, we obtain equations 
satisfied by the quantum components p^ of the density 
fields, 



dtPR + d x 



Ztt 



9_ 

2tt 



Pl 



d t pL -d x (v + —)p L + ^-Pr 

Z7T Z7T 



= J ■ 

= 0. 



(90) 



where the source term j(x,t) is defined by Eq. (|57p . The 
solution of Eq. (|9"0"]) determines the phases S^t) according 
to Eqs. (ffTlj) , (f6"3"| . Remarkably, Eq. (|6"5|) expresses the 
phase S v (t) affected by the electron-electron interaction, 
through the asymptotic behavior of p(x,t) in the non- 
interacting parts of the wire (regions I and III in Fig. 0| . 
The phases 8 v (t) determine the GFs via^ 



G%(r) 



where 



and 



. i A A R [S R (t)]A L [6 L (t)} 
'2ttu (liiAr^+T 



1= (1-K) 2 /2K, 



K = (1 + g/nv) 



-1/2 



(91) 



(92) 



(93) 



is the standard LL parameter in the interacting region. 

To explicitly evaluate 5 n (t) for the structure of Fig. |H 
it is convenient to rewrite Eqs. ([90]) as a second-order 
differential equation for the current 

J = v(p R -p L ), (94) 



(cu 2 + d x u 2 (x)d x )J(uj,x) =0, x t^O, (95) 
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where 



and 



u(x) = v(l + g(x)/irv 



,1/2 = 



K{x) 



(96) 



is a spatially dependent plasmon velocity. Reflection and 
transmission of plasmons on both boundaries is charac- 
terized by the coefficients r v , t„ (rj! + t"i = 1); here the 
subscripts rj refer to the boundaries between regions I/II 
and II/III. For simplicity, we assume them to be con- 
stant over a characteristic frequency ranged 4 u> ~ t _1 . 
The scattering matrices on the left and right boundaries 
have the form 



St 



TL t L 



(97) 



where the first component corresponds to the left mover 
and the second one to the right mover. 

Solution of Eq. (|9"5|) is quite straightforward. The 
boundary points x — ±L/2 and the observation point 
x = divide the x axis into four regions (I, II , 11+ , and 
III). In each of the regions the function J(uj,x) satisfies 
the homogeneous wave equation, with the velocity v (in 
regions I and III) or u (in regions II_ and 11+ ). The so- 
lution in each of the regions is thus a sum of two waves 
propagating left and right. As discussed after Eq. (fS"8")) . 
we need an advanced solution, which imposes the con- 
dition that in the leads (regions I and III) only incom- 
ing waves are present. There remain six coefficients that 
are fixed by the boundary conditions at the sample/lead 
boundaries [see Eq. (|9"T|) ] and by the matching condition 
at the observation point (x = 0). The latter condition is 
generated by the source term in Eq. (|9"U|) . 

Solving Eq. (f9"S")) and using Eq. (|M)) , we find the quan- 
tum density components p^. In accordance with Eq. (163|) 
the scattering phases S v (t) are determined by the behav- 
ior of p v in the asymptotic regions (x < —L/2 for pR and 
x > L/2 for p L ). We find 



p R {u,x) 



; | + K)tL e ikx + i ( k - K ) L / 2 
2V2Kv l-r R r L e~ 2iKL 



x ( 1 — r R re 



—inL\ 



1 



L 



1 



x<--; 



(98) 



_ (1 + K)t R e -ikx+Hk- K )L/2 
Pl[u,x) = == — tt—t 1 _ e 



x (r + r L e ikL ) , x > 



(99) 



Here we use the notations k — u)/v, n = u)/u, and r — 
(l-K)/(l+K). Substituting this in Eq. , we obtain 
8r,(t) in the form of a superposition of rectangular pulses, 



M*) = E 3v,nW T (t, t n ) . 



n=0 



where 



t n = (n + 1/2- l/2K)L/u 



(100) 



(101) 



5 v ,2m = 7rf_„r£V£(l + rjK)/VK 



S v ,2m+i = —irt- n r ri 



m+l m 



r™Jl-r]K)/VK . (102) 



For the "partial equilibrium" state (where n R (t) and 
n^ft) are of Fermi-Dirac form but with different tem- 
peratures and chemical potentials) the functional deter- 
minants are Gaussian functions of phases, reproducing 
earlier results of functional bosonization 28 . Indeed, us- 
ing Eq. (|A10j) . we find 



iA 



1 



2-ku (1 ± iAt) 1+ t 

l + Q 



■kT r t 



( ttT l t 



(103) 



\ sinh 7tTrt / \ sinh7r T^t / 
where the exponents 1 + a and (3 are given by the sums 



1 + a 



E 

oo 

E 

n=0 



OR.n 

~2~t7 

2tt 



(104) 



Substituting here the results (|102p for the phases 6 Vin , 
we obtain 



1 







Tr 



1 - K L K R 
Tr 



1 



7 



(1 + Kr) 



1 - TZ l TZr 



n L + ^-{i + n L ) 



(105) 



in agreement with Ref. |28| . One may check that, due to 
the sum rule 



+ /3 = 7, 



(106) 



at thermal equilibrium (Tr =Tl) the GFs are inde- 
pendent of plasmon transmission/reflection amplitudes. 

The phases 5^ (t) are shown in Fig. [5] for two limits of 
adiabatic {r n — 0) and sharp, 

r n = (l-K)/(l + K), 

boundaries. Let us stress that when we speak here about 
sharp boundaries, we mean that the extension of the con- 
tact regions is small compared to the characteristic plas- 
mon wave length ~ u/lu. It is assumed throughout the 
paper that the structure is always smooth on the scale 
of the electron wave length, so that no electron backscat- 
tering takes place. 

In physical terms (^(i) characterizes phase fluctuations 
in the leads that arrive at the measurement point x = 
during the time interval [0, r]. These fluctuations gov- 
ern the dephasing and the energy distribution of elec- 
trons encoded in the GFs G^(r). Up to inversion of 
time, one can think of S^it) as describing the fractional- 
ization of a phase pulse (electron-hole pair) injected into 
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FIG. 5: Phases 5 v (t) entering Eq. (fiJTj) for the GFs for sharp 
[a),b); r — (1 — K)/{l + K)] and adiabatic [c),d)] boundaries. 



the wire at point x during the time interval [0, r]. This 
is closely related to the physics of charge fractionaliza- 
tion discussed e arlieii2ii£2i£5iSSi£L££. At the first step, the 
pulse splits into two with relative amplitudes (1 + K)/2 
and (l — K)/2 carried by plasmons in opposite directions, 
cf. Refs. 31.66,671 As each of these pulses reaches the 
corresponding boundary, another fractionalization pro- 
cess takes place: a part of the pulse is transmitted into 
a lead, while the rest is reflected. The reflected pulse 
reaches the other boundary, is again fractionalized there, 
etc. Let us stress an important difference between bound- 
ary fractionalization of transmitted charg o 60 ! 68 and that 
of dipole pulses discussed here. While in the former case 
the boundaries can always be thought of as sharp (one 
is dealing with the small q limit), in the present prob- 
lem the way K(x) is turned on is crucially important for 
reflection coefficients r n at w ~ t _1 . 

For t <C L/u the coherence of plasmon scattering may 
be neglected and the result splits into a product 



(107) 



with each factor representing a contribution of a single 
phase pulse S n , n (t) = d v , n w T (t, 0). 

We now apply our general result (|9Tj) . (|107p to the 
"full non-equilibrium" case, when n^ie) have a double 
step form, Eq. (|69p. To obtain the exact form of the 
Green function G I) (r), one has to evaluate the Toeplitz 
determinants numerically. Here we restrict ourselves to 
the evaluation of the long-time asymptotics of (t) that 
can be found analytically employing Eq. (f66|) and governs 
the broadening of the split zero-bias-anomaly dip o 27 ' 28 . 
We focus on the adiabatic limit when the distribution 
function remains unchanged and the broadening is solely 
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FIG. 6: Dephasing rates 1/t^ r and ^/rf L as function of LL 
parameter K for the adiabatic case and double step distribu- 
tions with clr — cll — 1/3. 



due the non-equilibrium dephasing rat o 27 ' 28 , 1/V7. We 
obtain 



1/7^ = 1/7^ + 1/^ 



RL 



(108) 



where 1/t7 is the contribution to dephasing of the r\ 
fermions governed by the distribution of the rj' fermions. 
These dephasing rates are found to be 

2 vra + ^iT 



l/ r f7 = _£^l n 

1 * 2ir 



1 — 4a ?) (l — a v ) sin 



2VK 



(109) 

see Fig. |6] 

Two remarkable features of this result should be 
pointed out. First, let us compare our results with 
the results of RPA approximation. Consider the weak- 
interaction regime, 7 <C 1. We then obtain 

l/rf L ~njeV L a L (l-a L ), (110) 

l/rf R ~ ^( 7 2 /8)etWl - oh) • (Ill) 

This should be contrasted with RPA which predicts equal 
1/t^ l and 1/t£ r , see Ref.[27|. While l/r^ L agrees with 
the RPA result, X/t¥ r is parametrically smaller (sup- 
pressed by an extra factor of 7). The reason for this 
failure of RPA is clear from our analysis. For a weak 
interaction the contributions of R and L movers to Gr 
are given by the functional determinants A TtT (S v ) with 
phases (for adiabatic boundaries) Sl — (1 — K)ir and 
Sr ~ ir(l + K). While the contribution of the small phase 
Sl is captured correctly by RPA, a small-<5 expansion of 
In Ar t (Sr) fails for large Sr (apart from equilibrium and 
"partial equilibrium" where In A VT (S) oc 5 2 .) 

Another important observation is that for certain val- 
ues of the interaction parameter K (different for r\ = R 
and L) the dephasing rates I/t^ vanish. This implies 
that, for these values of K, the GF does not decay ex- 
ponentially in time, so that the power-law ZBA anomaly 
is not smeared. The absence of dephasing indicates that 
for these values of interaction the system reduces in some 
sense to a non-interacting model. As we are going to 
show, at these points the functional determinant can be 
calculated exactly. 



2. Refermionization 

The points of no-dephasing correspond to the value 
of the phase 5 (argument of the functional determinant) 
equal to 8 = 2im with an integer n. We will demonstrate 
that at these points the functional determinant A T (<5) 
can be calculated exactly by "refermionization". The 
case 5 = 2tt corresponds to the non-interacting (K = 
1) single-particle GF and has been already analyzed in 
Sec. IIII1 To study the case 5 = 4ir, we consider a two- 
fcrmion GF 



G 2 = (1^(1)^(2)^(3)^(4)) 



(112) 



We focus on the limit of merging points, t\ = t2 = 0, 
<3 = t± = t; Xi, X2, %3, Xi — ► x, which corresponds to 
simultaneous creation and annihilation of two fermions, 
and thus should generate A T (47r). For non- interacting 
electrons the GF G 2 can be readily calculated. Using 
Wick's theorem, we find 



G 2 = G(3, 1)G(4, 2) - G(4, 1)G(3, 2) 



(113) 



If the spatial points strictly coincide, x\ = X2 = X3 — 
x^ = x, the function G2 vanishes. A finite result is ob- 
tained after splitting the points by distances of the order 
of Fermi length, Sj ~ v/A. We thus find 



G 2 (r) 



-(- 
2 V 



where x% = x 
sponds to 



L - S 2 )(s 4 - S3) 

- d S2 ) 2 G(r, s\)G(t, s 2 ) |a 1= i s =o , (H4) 
Si. In the bosonic description this corre- 



= (JL\ 



2i0(r)-2i0(O) 



). (115) 



As was shown above, this correlation function can be 
evaluated, with the result expressed in terms of a func- 
tional determinant, 



G„(r) 



= {- 
\2ttv 



2 A t (4tt) 



(116) 



(At) 4 ' 
116|) . we express A r (47r) 



where we used t ^> A -1 . 

Comparing Eqs. (|114|) and 
through the free-electron GFs, 

A T (4^) = (2Tr) 2 (vTf(d Sl ~d S2 ) 2 G( Sl ,T)G( S2 ,T). (117) 

The numerical coefficient (2ir) 2 was restored by compari- 
son with equilibrium case. For a double-step distribution 
function, Eq. l[69|). we find from Eq. pT7|) : 



A r (47r) 



= e 2i(o-l)V 
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FIG. 7: Refermionization at a no-dephasing point. Upper 
plot: double-step distribution function 1 — n(e). Lower plot: 
function y"(e/V), where Gj 'M) = -*/(2^)(U/A) s i/(e/U) 
is the FES Green function (|77[1 proportional to the determi- 
nant A T (47r) with the argument 5 — 4ir being integer multiple 
of 2ir. The second derivative is plotted in order to emphasize 
singularities. The arrow at e = denotes a delta-function 
contribution to y" . 



We see that A T (47r) shows oscillations in r. At zero tem- 
perature there is no exponential damping. The absence 
of damping is a manifestation of the vanishing dephas- 
ing rate, see the discussion above. Another interesting 
property of the result (|118[) is the emergence of oscilla- 
tions with three frequencies: —2aeU, (1 — 2a)eU, and 
(2 — 2a)eV, implying three points of singular behavior in 
the energy space. Let us recall that the input double-step 
distribution had two such points: —aeV and (1 — a)eV. 
With increasing interaction strength the corresponding 
twofold singularity gets progressively more smeared [see 
Eq. 1(50])]. but then as 5 approaches 5 — 47r, a threefold 
singularity emerges at the new positions, see Fig. [7J 



(118) 



This procedure can be extended to a more general case 
of 6 = 2im. Indeed, the simultaneous creation and anni- 
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hilation of n non-interacting fermions is described by 

where we again imply a point splitting on a distance of 
the order of Fermi wave length, i.e. ~ v/A. (One can 
check that the relation resulting from this consideration 
does not depend on details of the point-splitting proce- 
dure). The function G n can be expressed in terms of 
single-particle GFs as follows: 

n 

G„(r) = C n s n ^l[(d St -d S] )G( Sl ,T) 

x G(s 2) T)...G(s n> T) U =0 . (120) 

Here C n are numerical coefficients of the order of unity. 
On the other hand, in the bosonic framework we have 



/a 

Gn[T) -\2^) (Ar)« 2 



(121) 



Demanding the equivalence of Eqs. (|120[) and (|121[) . we 
establish the identity 

n 

A T (2nn) = C n (vrf Y[(d Si ~ 9 'i ) 

x G(s 1 ,T)...G(s n ,T)) | s=0 , (122) 

expressing the functional determinant A T (<5) through free 
fermionic GFs G(r) for S = 27m. The numerical coeffi- 
cients C n can be restored by comparison with the known 
result for A T (27m) at equilibrium. The explicit form of 
A r ((5 = 2nn) can be readily found by substituting in 
Eq. p22p an explicit expression for the GF for a given 
distribution function. 



3. Tunneling into non-interacting regions 

Next we discuss the tunneling spectroscopy for the 
non-interacting parts of the wire. Let us focus on the 
right-moving electrons; the analysis of left-moving ones 
can be done in the same way. For x\,x 2 < —L/2 (re- 
gion I in Fig. the GF is the one of free fermions, as 
the right-moving particles emerging from the left reser- 
voir are not yet aware of the interacting region they 
are about to enter. The situation is less trivial for 
xi,X2 > L/2 (region III). Indeed, while the strength of 
the interaction in this region is zero, right-moving elec- 
trons there have passed through the interacting part of 
the wire, which modifies their Green function. We will 
show below that the GFs G^(x\ 1 1\; x 2 , t 2 ) in the non- 
interacting region satisfy Galilean invariance: they de- 
pend on {x\ — vti) — (x 2 — vt 2 ) only. For this reason, it 
is sufficient to consider x\ — x 2 to obtain the full infor- 
mation about the GF. 



The evaluation of the GF is performed in the same way 
as in the interacting region, yielding the result (ji?Tj) . The 
phases 8r,{t) are now given by 



**(*) = J2S R , n w T (t,x 1 /v + 2t n ), (123) 
x\ — L N 



n=0 

S L (t) = 2nr R w T ( t, 



( x L \ 

+ V S L , n W T I t, — + - + 2tn J , (124) 

with the following amplitudes of rectangular pulses: 



S Rin = 2nt L t R (r L r R ) n , 
Sl,u = -2n{r L r R ) n r L t R . 



(125) 



In the case of smooth boundaries only one pulse is cre- 
ated, 5 R ,o = 2ir, reproducing the free fermion GF. Thus, 
in the adiabatic case, the interaction has no influence 
on GFs in the non-interacting parts of the wire, as ex- 
pected. If the transition between non-interacting and 
interacting parts of the wire is not smooth, plasmon scat- 
tering takes place. This process leads to a redistribution 
of electrons over energies^ and thus affects GFs in the 
non-interacting region. 



C. Green functions at different points and 
Aharonov-Bohm interferometry 

So far we have discussed GFs at coinciding spatial 
points, having in mind tunneling spectroscopy experi- 
ments. We now consider GFs at different spatial coordi- 
nates. Such GFs are relevant to various physical quan- 
tities, in particular, in the context of Aharonov-Bohm 
interferometry. The similar problem in the context of 
chiral edge state has been considered in Refs. [69lf70l Let 
us consider a four-terminal setup formed by two quantum 
wires coupled by tunneling at two points, as schemati- 
cally shown in Fig. [SJ Each one of the quantum wires is 
assumed to be a LL conductor connected to two non- 
interacting electrodes with arbitrary (in general, non- 
equilibrium) distribution functions, as shown in Fig|31 
We are interested in the Aharonov-Bohm effect, i.e. the 
dependence on the magnetic flux $ of the electric cur- 
rent flowing from wire 1 into wire 2. Consider the situa- 
tion where the tunnel coupling between the wires 1 and 
2 is weak. We also assume that both arms of the AB- 
interferometer have equal length d and tunneling occurs 
at points located inside the interacting part of the wire. 
The flux dependent part of the electric current is given 
by 



G 2 > (d,t)Gf(-d,-t)]+h.c. 



(126) 
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FIG. 8: Aharonov-Bohm setup in a four terminal geometry, 
with tunnel coupling (dashed lines) at two points. Both in- 
terferometer arms have length d. 



where the subscripts 1 and 2 label the wire and t\2 is the 
tunneling matrix element between the wires. Separating 
the GF into left and right moving part, one gets 

/oo 
dte-^[G< 7j (d 7 t)G> v (-d } -t) 
-oo 



G>Jd,t)G<J-d,-t)}+h.c. 



(127) 



where we have neglected terms that oscillate fast with 
the interferometer size d. 

To analyze the GF G^ (x\, X2, t) between two different 
points of a wire, we proceed in the same way as in the case 
x\ = X2 above. Integration over the classical component 
of the density field leads to equations of motion for its 
quantum component (we choose ?/ = R for definiteness) , 

(— ioj + vd x )p R<u +d x (J^-Pu} = j{u,x;xi,x 2 ,T) , 
{iuj + vd x )p LtU) + d x (£-p u ) = , (128) 



where we have used the (uj,x) representation; p u — 
P~r,l> + P~l,uj- Equations ()128|) differ from the earlier 
Eq. (|90p only by the source term, which now reads 



j(uj,x;x 1 ,x 2 ,T) 



1 

72 



6(x-x 1 )e iwT -6(x-x 2 ) 



(129) 



Solving Eq. (fl28|) . we find for x < -L/2 

(1 + K)t L e ikx +i( k - K ) L / 2 



Pr{u,x) 



-2inL„ 



e -inx 2 _ e i{wT-KXx) _ rRre il y x ^^ L ) K 



r R re 



jr-\-i(xi —L)k 



(130) 



Similarly, we find for x > L/2 

(1 + K)t R e-ikx+iik-^L/ 2 



2y/2Kv l-e- 2iKL r R r L 



iur — {xi-\-L)K 



r L e 



—i(x-2+L)n 



(131) 



Employing Eqs. ((Ml), QUI and ([I3T]) . we obtain the 
following result for the GF: 



G^{x 1 ,x 2 ,t) = - 



1 



27ra(±iA)T 
A R [5 R (t)]A L [S L (t)} 



Xl-X 2 -r- J_\l + " 

u * A' 



(132) 



It is interesting to note that for spatially separated points 
the scaling of GF with time (and consequently with en- 
ergy) is affected by plasmon scattering at the boundaries 
between wire and the leads. Surprisingly, even at equi- 
librium the GF inside interacting region is affected by 
the way interaction is turned on. For coinciding spatial 
points the universal LL exponents, characteristic of an 
infinite wire, are restored due to the sum rule (I106|) . 

In the long wire limit the functional determinant splits, 
as before, into a product 

oo 
n=0 

(133) 

Here S v , n are given by Eq. (|102[) . The calculation of 
is performed in a similar way, yielding 

oo 

A L [S L (t)] ~ Yi ,th »i-'2^r,.2n)A £| . ^^{8lm+i)- 

n=0 

(134) 

We see that in the case of a GF at different spatial 
points the time argument of A ?7iT (5 r) ) (determining the 
duration of the pulses) is replaced as compared to the 
case of x\ = X2 by 



r T V- 



X\ - X 2 



(135) 



with the — (+) sign corresponding to even (respectively, 
odd) pulses. It is easy to understand the reason for these 
alternating signs. The even pulses are those that experi- 
ence an even number of reflections, thus preserving their 
chirality, while the odd pulses experience an odd number 
of reflections and thus invert their chirality. We note that 
in the case when both points are located in one of the non- 
interacting regions (x±,X2 > L/2 or X\,X2 < —L/2), the 
same consideration leads to an analogous replacement of 
the time argument but with the bare velocity v, 



Xi - x 2 



(136) 
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in the phases S v (t), Eqs. (| 1 23|) and (|124[) , entering 
A^ T (5 I) ) and, correspondingly, in GFs. Since there is no 
fractionalization at the tunneling processes into a non- 
interacting region, only half of the pulses survives and 
no sign alternation arises. 

Of particular interest is the value of GF at the 
interaction-renormalized "light cone", x\ — X2 — ±ut. 
The value of the GF at these points determines the in- 
tegral for the interference current in Eq. (|127[) . see also 
Refs. 12211311 . Let us consider for simplicity the case of adi- 
abatic barrier (r^ = r R = 0) when each of the products 
(|133p . (|134[) reduces to the first factor. Compared to the 
limit of coinciding spatial points, the duration of pulse in 
the functional determinant has changed. The contribu- 
tion associated with (x\ ~ X2 — r/ut) leads to a doubling 
of pulse duration in , while A,, has disappeared al- 
together. As we see now, the dephasing rate governing 
the exponential damping of the GF at x\ — xi = tj'ut 
is 



lf W) 



(137) 



where 1 /t^ are the partial dephasing rates for the tun- 
neling spectroscopy problem (coinciding spatial points), 
as introduced in Sec. IV B II The dephasing rates (|137p 
manifest themselves in the interferometry measurements 
by inducing an exponential damping of the corresponding 
contributions to the Aharonov-Bohm oscillations (thus 
the superscript "AB" ) . In the limit of large interferom- 
eter size d, the contribution with the lowest dephasing 
rate will dominate, 



AB 



min 2/T7" 



(138) 



For double step distributions, the dephasing rates 
([T5T)) are given (up to a factor of 2) by Eq. (fTUi?)) . With 
increasing interaction the dephasing rate 1/t£ b begins 
to oscillate as a function of interaction parameter K, as 
illustrated (for the case of adiabatic contacts with leads) 
in Fig. [5] This leads to a remarkable prediction: the visi- 
bility of Aharonov-Bohm oscillations should be a strongly 
oscillating function of the interaction strength. 



D. Spinful Luttinger liquid 

We now consider the problem of tunneling spec- 
troscopy for spinful electrons. The analysis is a straight- 
forward extension of the spinless case, analyzed in 
Sec. lVBl We begin with a fermionic Hamiltonian, which, 
in the spinful case, is given by 

H = H a + H ce , (139) 

H Q = -ivY,{^R,*d^R,<r-A,ML,*) ,(140) 



(141) 



where the index a =T,| labels the spin projection. We 
now switch to a Lagrangian description. To construct the 
free part of the action on the Keldysh contour we repeat 
the steps described in detail in Sec. IIIII and find 



So = Pn.^K W-tk-Zfcwr], ( 142 ) 

where 

X v , a = p v . a . (143) 
The interacting part of the action reads 



j dxg{x)p v ^p 7] ,y 



(144) 



To describe the tunneling spectroscopy measurements, 
we need to find the single-particle GFs 

Gl v A x > f ) = -WvA*> *)V4,a(0, 0)> • (145) 

The fermionic operators are expressed in terms of bosonic 
fields (which now also carry the spin label) in the usual 
way, 



(— 

\2irv 



1/2 



(146) 



Substituting Eq. (|146p into Eq. (|145|) , representing the 
GF as a bosonic functional integral with the action 
Sq + S ee , and performing the integration over the classi- 
cal component of the density field, we find the equation 
of motion satisfied by the quantum components of the 
field, 

(d t + vd x )p R + d x ^-(p R + p L ) = j(x, t) , 
(-d t + vd x )p L +d x ^-(p R + p L ) = 0, (147) 

and 



(d t + vd x )s R = j(x,t) , 
(-dt + vd x )s L = 0. 



(148) 



Here we have passed to new variables that describe the 
spin and charge sectors of excitations, 

PR = PRA + PB„i , Pl — PL,] + Pl, i 

sr = Pr,] - p R ,i , s L = p L ,T - PL, i ■ (149) 

As one sees, the equations for the charge and spin de- 
grees of freedom are decoupled, which is a manifestation 
of spin-charge separation. The spin-density component 
obeys the same equation as the density of free fermions, 
Eq. Therefore, the spin sector is characterized by 

a LL parameter K s = 1. As follows from Eq. (|148p . the 
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spin component s v propagates through the wire without 
any reflection. 

To find the charge component, we define the charge 
density current 



J = v(PR - PL) ■ 



(150) 



In terms of J, equations (|147p are reduced to a second 
order differential equation, 



(uj 2 + d x u 2 c d x )J = for x ^ 0, (151) 



where 



ui(x) = v A /Kt(x), 



K n = 1 



2.9 



-1/2 



(152) 



Equation (|151[) coincides with Eq. (|9"B")) up to a different 
definition of the LL parameter. The interaction parame- 
ter 7 = (1— K C ) 2 /2K C and the transmission and reflection 
amplitudes are determined as for spinless fermions, with 
the replacement K — > A' c . 

The resulting expression for the Green function of spin- 
ful fermions within the non-equilibrium bosonization ap- 
proach reads 



iA 



ruA^[<v(*)] 



2tt 



(1 ± iAr) 1 



H7/2 



(153) 



Here we have assumed for generality that distribution 
functions of spin-up and spin-down particles may be dif- 
ferent. Therefore, the distribution function is labeled by 
two indices (chirality and spin projection); these indices 
are inherited by the functional determinant. The time- 
dependent phases of the spinful fermions 5, ha (t) are ex- 
pressed in terms of the scattering phases <L (t) of spinless 
fermions, Eq. (jlOOp , in the following way: 



*l,t(*) = *£,;(*) = 5*i(*). ( 154 ) 

**,t(*) = 5 (<**(*) > (155) 
6R,l(t) = \(6R(t)-#k(t)) , (156) 

where 5 R (t) corresponds to non-interacting fermions and 
consists of a single pulse with an amplitude 2ir, S R (t) — 

2TTW T (t,0). 

We conclude that the inclusion of spin changes the 
scattering phases in an essential way. This is most im- 
portantly seen when considering the first pulse propa- 
gating to the right. Let us assume that there is no 
reflection at the boundaries with non- interacting leads. 
The corresponding scattering phases are each a super- 
position of the spin and charge modes, see Eqs. (|155[) . 
(|156p . Since the velocities of these modes are different 
(v and it, respectively), then for sufficiently long wires, 
L(u _1 — u _1 )/t ^> 1, the first pulse splits into a charge 



and a spin parts. For a short wire, the spin pulse and the 
first charge pulse overlap. In this case one has to deal 
with the general formula (|153[) and with time-dependent 
phase containing both, spin and, charge contributions. 
Hence, if the wire is sufficiently short (or, in other words, 
for a given length of the wire the interaction is sufficiently 
weak) the spin-charge separation does not have enough 
time to develop. For sufficiently long wires, spin-charge 
separation does take place, in which case the respec- 
tive determinants can be written as products of spin and 
charge contributions. This decomposition is not valid for 
short wires (or, for a given length of the wire, for suffi- 
ciently weak interaction). Note that at equilibrium there 
is significant simplification. The GFs depend only on 
the sum of the scattering phases squared, see Eq. (|103p . 
Due to the sum rule (|106|) this combination remains un- 
changed, regardless of whether spin and charge pulses 
overlap or not. Thus, at equilibrium one can always 
think about these two modes separately and additively. 
Out of equilibrium, the dependence of the GF on scatter- 
ing phases is more subtle, and the results for overlapping 
and separated pulses are different. Therefore, the spin- 
charge separation occurs in this case only for sufficiently 
long wires. Focusing on this regime, we find 



A fi , T [fe, T ] = A fl , r , T (7r) Yl Afl,r, T 

n=0 
oo 

AiU^flJ = Afl, Tli (-7r) JJ Ah,t 

OO 

A L . a [SL.a} = n^Lr 

n=0 



Sr 



n=0 



, (157) 



,(158) 



(159) 



The first factor in each of Eqs. (|157p and (|158p originate 
from the spin mode yielding the phase tt, i.e. a half of 
the free-fermion phase value. The scattering phases of 
other pulses (originating from the charge mode) have a 
half of their values for spinless electrons. 

Let us analyze this result. Consider the case of smooth 
(adiabatic) contacts with leads, so that only one pulse 
passes in each direction (all <L „ with n > 1 are zero). 
For the case of partial equilibrium, the determinants can 
be evaluated explicitly, yielding 



n A R As R At)} = ( 



nT R T 



V sinh 7tTrt 



l+a/2 



I1^M0] = (-^t) , (160) 



\sinh ttTlt 



where a,f3 are given by Eq. (|105p . For a double-step 
distribution function, a scmiclassical limit of the determi- 
nants (|157p . (|158p . (|159p can be readily evaluated. In the 
large-time limit the behavior of the r.h.s. of Eqs. (|157p . 
(I158p . (|159p is exponential, yielding the partial decay 



rates 



^RR 



eV 

7T 



In (1 - 4a(l - a)) 



n=0 

eV 

7T 



5^ In ( 1 - 4a(l - a) sin' 



n=0 



In [ 1 - 4a(l - a) shr 



2 ^i,n 



(161) 



and the total rate T R = T RR + T HL . Let us stress an 
important difference with the spinless case. There, for 
smooth boundaries, the distribution function was not 
affected, and the decay rate (inducing the smearing of 
singularities in tunneling spectroscopy) was solely due 
to dephasing. In the spinful case the situation is dif- 
ferent: independently of the shape of the boundary the 
spin-charge separation affects the distribution function of 
electrons. Indeed, imagine that we perform the tunnel- 
ing spectroscopy of right-movers in the right lead (non- 
interacting region III of Fig. [4]) for the case of adiabatic 
boundaries. Then the phases are $L,n = $R,n=to = and 
<5r.o = 2-7T. In the spinless case this implied that the dis- 
tribution function remained unchanged. This is not so in 
the spinful situation, however: according to Eqs. (|157l) . 
(|158p we get now a product of four determinants with 
arguments ±7r, 



[A fliT:T W] 2 A^ ra (-7r)A fi , Tii W 



(162) 



implying that the distribution function has changed. 
This effect remains finite even in the limit of vanish- 
ing interaction (K — > 1) as long as the long- wire con- 
dition, LTiuT 1 — v^ 1 ) ^> 1, is satisfied. Returning to the 
spectroscopy of the interacting region, we conclude that 
both effects — dephasing and change of the distribution 
function — are necessarily present in the spinful LL case 
and cannot be easily "disentangled" . The decay rates T 
presented in Eq. (|161[) and in Fig. [!|] yield the combined 

effect of interaction on the GF G\ (t) and determine the 
smearing of tunneling spectroscopy singularities in the 
energy space. 

Finally we discuss the extension of Eq. (|153[) to the 
case of GFs at different points. In this case we find 
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FIG. 9: Decay rates r RR and T RL (governing smearing of 
tunneling spectroscopy singularities) as functions of LL pa- 
rameter K c for spinful fermions. The adiabatic coupling to 
leads and the double-step distributions with an = cll — 1/3 
are assumed. 

where x = x\ — X2 and r ? '' = r + r/(— l) n+1 x/u. As for the 
spinless case, Eq. (|132[) . the scaling of GFs with spatially 
separated points is affected by plasmon scattering at the 
boundaries between the interacting regions and the leads 
even at equilibrium. 

Before concluding this section, we point out a connec- 
tion between the spinful LL and the problem of the in- 
teger quantum Hall edge with two edge channels (corre- 
sponding to two Landau levels below the Fermi energy 
in the bulk). Non-equilibrium properties and quantum 
coherence of such a system are currently attracting large 
interest, in particular, in connection with experiments 
on quantum Hall Mach-Zehnder interferometers^. In 
the quantum Hall setup the edge channel index plays 
a role of spin. The main difference is that, under con- 
ventional circumstances (if one does not make special 
efforts to couple counterpropagating edge modes), the 
quantum Hall system is chiral: there are, say, only right- 
moving modes and no left-movers. This leads to a num- 
ber of essential simplifications: (i) the tunneling density 
of states becomes trivial (no ZBA); (ii) the charge frac- 
tionalization is absent (since plasmons can move only in 
one direction). What remains is the charge-spin sepa- 
ration. This implies that following simplifications with 
functional determinants (|157[) . (|158p . (| 1 50|) the product 
of which determined the tunneling spectroscopy Green 
function G^^(t) within our analysis. First, the left- 
mover determinants (|159[) are now absent. Second, out of 
the set of phases 5n^ n only the n — phase remains, be- 
ing equal to its free-fermion value, Sr.q — 2n. Therefore, 
the product of determinants takes the form (| 1 62[) (that 
has appeared above in the context of a GF in the non- 
interacting region of a non-chiral spinful wire with adia- 
batic contacts). The analogous statement holds for the 
Green function with different spatial points, Eq. (| 1 63[) . 
Specifically, for the case of the two-channel chiral setup 
(relevant in the quantum Hall Mach-Zehnder interferom- 
etry context^) the last fraction (having L determinants 
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in the numerator) in Eq. (| 1 63|) disappears, while in the 
preceding-to-it fraction one should keep only n — fac- 
tor and set £r i0 = 2tt. An equivalent result was obtained 
by a different method in the recent work^. 



VI. SUMMARY 

Let us summarize the main results of this work, fol- 
lowing the flow of our presentation in the paper. 

1. We have developed a nonequilibrium bosonization 
approach and derived a bosonic theory describing 
the LL of interacting ID electrons out of equilib- 
rium. The theory is characterized by an action de- 
pending on density fields defined on the Keldysh 
time contour. In contrast to the equilibrium case, 
this theory is not Gaussian, which is a manifes- 
tation of the fact that the density matrix is non- 
diagonal in the bosonic Fock space. We have used 
this theory to calculate the electronic GFs govern- 
ing various observables. 

2. We have first calculated the GF of non-interacting 
fermions from our non-equlibrium bosonization ap- 
proach. The GF is expressed in terms of a func- 
tional determinant of the Fredholm (more specifi- 
cally, Toeplitz) type similar to those that have ear- 
lier appeared in the context of counting statistics. 
The key difference is that in the case of counting 
statistics the determinant is non-analytic and 2tt- 
periodic in the counting field (which reflects charge 
quantization), while in our theory the determinant 
should be understood as an analytically continued 
function. We have found that the free-fermion GF 
is described by the determinant exactly at the point 
2ir, which is related to the fact that in the bosonic 
theory a fermion is represented by a 27r-soliton. 

3. We have next generalized the GF calculation to the 
problem of non-equilibrium Fermi edge singularity 
describing excitation of an electron into the conduc- 
tion band within the process of photon absorption, 
accompanied by creation of a core hole. The result 
is obtained in terms of the same functional determi- 
nant as in the free-fermion case but the argument is 
now shifted from 2tt by twice the scattering phase 
on the core hole. 

4. We have then applied our formalism to the problem 
of interacting ID fermions. We have considered a 
model of a LL wire coupled to non-interacting ID 
leads, with the interaction strength "turned on" in 
specified fashion at the boundary between the wire 
and each of the leads. We have shown that the elec- 
tron GFs — which describe tunneling spectroscopy 
measurements — are again expressed in terms of 
Fredholm determinants. The phases S^it) enter- 
ing the expressions for the corresponding operators 



have a physical interpretation in terms of fraction- 
alization processes taking place during the tunnel- 
ing event, near the boundaries. If the characteris- 
tic energy scales for the tunneling spectroscopy are 
large compared to the inverse flight time through 
the LL wire (Thouless energy) — which means that 
we are considering the truly ID (rather than 0D) 
regime — the functions S^it) represent a sequence 
of rectangular pulses separated by large intervals. 
As a result, the Fredholm determinant splits into a 
product of Toeplitz determinants of the same type 
as in the cases of non-interacting fermions and the 
Fermi edge singularity. 

5. We have analyzed the long-time asymptotics of the 
determinant which yields the dephasing rate con- 
trolling the smearing of LL tunneling singularities 
(zero-bias anomaly). The dephasing rate for the 
GF of electrons with rj (±1) chirality is a sum 

of two terms Y] n'—±i ^-/ T 2 originating from func- 
tional determinants which depend on the distribu- 
tion function of left- (rj' = —I) and right- (rj' = 1) 
moving electrons, respectively. For the case of 
double-step distributions, there are two important 
findings: 

(i) At weak interaction, comparing our exact results 
with those of the RPA, we find that while 1/tJ' ~ v is 
correctly obtained (to leading order) within RPA, 
the RPA result for I/t^' 77 is parametrically wrong. 
This demonstrates that even for a weak interaction 
a naive perturbative expansion (leading to RPA) 
may be parametrically incorrect in LL out of equi- 
librium. 

(ii) Both 1/tJ' ±v are oscillatory functions of the in- 
teraction strength (or, cquivalently, LL parameter 
K). Furthermore, each of them vanishes at certain 
values of K . At these values the "counting phase" 
for the corresponding determinant becomes an inte- 
ger multiple of 2n. We have calculated the determi- 
nants at these no-dephasing points by a refermion- 
ization procedure. 

6. We have generalized the above results to the case of 
a GF with two different spatial arguments. When 
considering the value of the GF G v at its main 
peak, x\ — X2 = rjut, the dephasing rate is 2/t^' - ' 7 , 
while I/tJ'* 1 does not contribute (and thus RPA 
is restored for weak interaction). The situation is 
reversed for the value of the G v at the other peak, 
x\ — X2 = —rjut, where the dephasing rate is 2/t^' 11 . 
Such GFs (with x\ — X2 — irjut) enter the expres- 
sion for the interference contribution to current in 
an Aharonov-Bohm interferometer formed by two 
LLs coupled by tunneling at two points. Our re- 
sults imply that the dephasing rate in such a non- 
equilibrium LL interferometer (and thus the visibil- 
ity of Aharonov-Bohm oscillations) is an oscillatory 
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function of the interaction strength. 

7. We have also considered the case of a spinful LL. 
The general structure of the results for the GFs 
is similar, with the key difference being that now 
we encounter products of determinants with phase 
arguments corresponding to the spin and charge 
sectors. This is a manifestation of the spin-charge 
separation. One important consequence is that the 
temporal decay rate of the Green function (and 
thus the smearing of singularities in the tunneling 
spectroscopy) remains finite in the limit of vanish- 
ing interaction strength, assuming the limit of the 
large system size is taken first. With increasing 
interaction strength, the dephasing rate oscillates, 
similarly to the case of spinless fermions. 

The non-equilibrium bosonization formalism devel- 
oped in this work has a variety of further applications. 
They include, in particular, counting statistics of charge 
transfer in an interacting ID system away from equi- 
librium, analysis of many-body entanglement, quantum 
wires with several channels, etc. Generalizations or mod- 
ifications of our formalism should be useful for a number 
of further prospective research directions, such as sys- 
tems of cold atoms and fractional quantum Hall edges 
away from equilibrium. 



VII. 
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Explicitly calculating the correlations functions of the 
bosonic fields one finds 



doj 



of A 



(1 — cos ujt) coth — ± i sinwr 



(A3) 



Next we calculate the integrals that appear on the r.h.s. 
of Eq. (|A3[) . The integral with sin ujt yields 



e U / A sin ujt — arctan At . 



(A4) 



The remaining integral can be split into two parts: 

50 UJ 

>o 

and 
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(A6) 



in the second one we have dropped the convergence fac- 
tor, e~ u < A , which is justified in view of T <C A. Employ- 
ing Eqs. (|A4| . l[A5j) . and (|A6]) . one gets 



J*(r) 



In 



/ ttTt 
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\ sinh -kTt 1 ± it A 
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Substituting Eq. (|A7|) into Eq. (jAip . one recovers the 
result for the GF of free fermions, Eq. (|54[) . 
We proceed with GFs for FES at equilibrium, 



G*(t) 



iA 

T 2^ 6XP 



{(1- 5oM 2 J^(t)} . (A8) 



Next we relate the GF and the functional determinant 
A T (<5). At equilibrium the latter can be evaluated as 
follows: 



APPENDIX A: EQUILIBRIUM: GREEN 
FUNCTIONS G,| , G% VIA BOSONIZATION AND 
FREDHOLM DETERMINANTS A T 

GFs of free fermions at equilibrium can be readily 
found. Since at equilibrium the bosonic action is Gaus- 
sian, the functional integration over bosonic fields is 
straightforward. The fermionic GFs is thus expressed 
as 



G < (r)=T^e 
Zirv 



(t) 



(Al) 



in terms of the bosonic correlation functions 

J> (r) = ~{T K [i>+, v (0, 0) - ^-,1,(0, r)] 2 ) , 
J K (r) = ~(T K [^ !V (0,0) - 0+,,(O,r)] 2 ) . (A2) 
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Using Eq. (|A4|) . we find 
A T (<5) 



sinh tttT J 



{1 + t 2 A 2 )^ & / 2 ^ 2 



(A9) 



(A10) 



Comparing Eq. (|A10|) with Eq. (|a7|) . we establish the 
exact relation (including the proportionality factor) be- 
tween the GF and the functional determinant, 



G<(T) ^ T 2^ \T±IAt 



i(l-So/n) 2 



A T (2ir- 26 ). 

(All) 
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We notice that the determinant A r (<5) as given by 
Eq. (|A10[) is a product of the temperature dependent 
and independent parts. It is convenient to normalize the 
result by its zero temperature value, 



(1 + t 2 A 2 )^ s / 2 * 
We thus present Eq. (|A10|) in the form 

A T (S) 



( 1 + t 2 A 2^(S/2tt)2 ' 

where A T (S) is the normalized determinant, 



(A12) 



(A13) 



M*) = (-7 



tttT 



\ (5/271 



V sinh 7rrT . 



(A14) 



By construction, A T (5) = 1 for T = 0. It turns out 
to be more convenient to deal with the normalized de- 
terminant, since all ultraviolet divergences (A-dependent 
factor) are excluded from this quantity. 



APPENDIX B: HIGH-ORDER VERTICES FOR 
ID FERMIONS: DIAGRAMMATICS 

In this Appendix we briefly sketch an explicit calcu- 
lation of third-order fermionic vertices by means of di- 
agramatic fermionic approach. Consider the third-order 
vertex shown in Figj2]3, 

5 3 ^(1,2,3) = (7^(1)^(2)^(3)) 

= (^4,7,^^,^2,77^3,^3,77) • (Bl) 

Here the index i — 1,2,3 includes the corresponding spa- 
tial coordinate (a;,), time (U) and Keldysh index Si that 



labels upper and lower branches, 
we find 



Using Wick theorem, 



-iS„(l,2,3) = G„(l, 3)^(3, 2)G,(2,1) 
+ G,(1,2)G,(2,3)G,(3,1). 

We choose first the following combination of Keldysh in- 
dices: si = +,s 2 = +,S3 = — . Using Eq. (fl6|) and 
passing into energy-momentum representation, we find 

£3,77(^1, qi, +; w 2 , q2, +; w 3 , 93, -) 

= --^tt) 2 ^)^) 
de 



x [1 - n ?) (e + lji) - n n (e + uj 2 )} . 



(B2) 



Here A t = u t - -qvqi and ui 3 = -u>i - w 2 , 93 = -?i - <72- 
Therefore, the third-order correlation function is re- 
stricted to the light-cone with respect to all its coordi- 
nates. At equilibrium the integration over energy yields 
zero, making the correlation function vanish. On the 
other hand, in the non-equilibrium situation the result 
is in general non-zero. Repeating the calculations for all 
possible choices of Keldysh indices s\, Sa, S3, we find that 
all third-order vertices are equal and given by Eq. (|B2[) . 
When transformed from s = ± to quantum and classical 
components, this means that the result is non-zero only 
when all indices are classical. Thus, the explicit diagram- 
matic calculations of the third-order vertex confirms that 
(i) it is restricted to the light-cone, and (ii) only corre- 
lations of classical fields are non-zero. We have checked 
this also for vertices of fourth order. These results are 
in full agreement with the general treatment (valid for 
vertices of all orders) performed in Sec. IIII1 
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